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1. Introduction 

Matching factors are needed to convert matrix elements calculated on the lattice into physical 
observables in continuum renormalization schemes. We focus here on matching factors for the 
bilinear operators that we are using as part of an ongoing numerical study using MILC configura- 
tions. The valence fermions of which these operators are composed are HYP-smeared ^ improved 
staggered fermions, which thus differ from the sea quarks, which use the "asqtad" staggered ac- 
tion [p. We calculate the matching factors (which are, in general, matrices) using perturbation 
theory at one-loop order. What is new compared to previous calculations using HYP-smeared stag- 
gered fermions [|3|, 0, ^, ^ ^ ^] is the use of the Symanzik-improved gluon action rather than the 
Wilson plaquette action. This generalization is necessary because the MILC configurations use an 
improved gluon action. 

This work is part of a larger project, whose ultimate aim is to provide matching factors needed 
for the calculation of both quark masses [Q] and Bk [10, 11, 12, |T^. The results presented here 



allow us to study the relative impact of improving quark and gluon actions. 

2. Actions and Bilinear Operators 

Since sea-quarks enter first at the two-loop level (through vacuum polarization diagrams), we 
do not give the details of the sea-quark action. 

For the valence quarks, we use the HYP-smeared staggered fermion action. This has the same 
form as the unimproved staggered fermion action, 

= ^L%(")^(")(^M(«)^(« + A)-'t^/I(«-A)z(«-/t)) +mY,X{n)X{n), (2.1) 

[with n G a lattice coordinate and T]^ {n) = (— l)"l+"2^ i""fi-i] but with the original thin links 
replaced with HYP-smeared links V^. For the details of HYP-smearing, see Ref. [[I]]. This action 
has a number of important properties: 



1. It substantially reduces the breaking of taste symmetry [ [141 , p3| , g|; 

2. It significantly reduces one loop corrections [^ ; 

3. It performs tadpole improvement automatically [|l6||; 

4. It reduces scaling violations efficiently [|T7|]. 

Given also its simplicity, it is an attractive choice for an improved staggered action. 
The HYP-smeared links can be expressed in terms of blocked gauge fields : 

{x) = exp {iB^ ix+fl/2)) . (2.2) 

The blocked gauge fields B^ can be written as a power series in the original gauge fields Ay, 

oo 

5m = L^?' (2-3) 

n=l 
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in] 

where contains all terms with n powers of the gauge fields Ay. In the present one-loop cal- 
culation, we need only the linear term, B^'. One might have expected the quadratic term B^^'^ 
to enter as well, but its contribution vanishes due to the fact that HYP-smearing includes SU (3) 
projections p^. 

The linear term can be written in terms of a blocking kernel h^y, which it is convenient to 
define in momentum space: 

(2.4) 



V 

All information about the HYP-smearing, including the smearing parameters, is contained in h^y. 
Following Ref. |^, we decompose the kernel into diagonal and off-diagonal parts 

h^v{k) = 5^yD^{k) + (1 - 5^y)G^y{k) , (2.5) 

With smearing coefficients chosen to remove &{a^) taste symmetry breaking coupling at tree level 
(«! = 0.75, a2 = 0.6 and a-i = 0.3 in the notation of Ref. [|l|]), the diagonal part is 



D,{k) = l-l^sl+ £ 



with = sin(^^/2), while the off-diagonal part is 



v<p 



-2 -2 



-2 -2 -2 

- SySpS^ . 



(2.6) 



-2 -2 

sis: 



■ + ■ 



(2.7) 



For unimproved staggered fermions the blocking kernel simply reduces to h^y = d^y. 
We now turn to the Symanzik-improved gluon action, which can be written as 



Cpi^ReTr(l - U^i) + CnJ^^eTril - +Cpg£ReTr(l - t/pg) 

pi rt pg 



(2.8) 



Here, pi, rt, and pg represent plaquette, rectangle and parallelogram, respectively. The coefficients 
c,, where /=pl, rt, or pg, should be chosen so as to improve the scaling behavior. In the present 
calculation, we use the tree-level improved coefficients 

5 1 



3' 



Crt 



12' 



and 



'Pg 



0, 



(2.9) 



which were determined by Luscher and Weisz QlSj [19p. The MILC collaboration actually use the 
1-loop improved values for the coefficients, but this leads to changes which would enter only into 
a two-loop calculation. Note that the Wilson gauge action corresponds to Cpi = 1 , Cn = Cpg = 0. 

For the lattice bilinears, we use the operators of Ref. [pO|], which reside on 2^ hypercube. The 
operator with spin 5 and taste F can be written as 

1 



[S X F]{y) = —'£[xb(y+A) {Ys(^^f)ab Xc{y + B)] r''{y+A,y + B) , (2.10) 



A,B 



where y is the coordinate of hypercube, and A and B denote the corners of hypercubes. Gauge 
invariance is maintained by the inclusion of '^'''^{y + A,y + B), which is constructed by averaging 
over the shortest paths connecting y+A and y+B, with each path constructed from the products of 
HYP-smeared links . In this way the operators are improved in the same fashion as the action. 
This also ensures that the currents [V x S] and [A x P] are conserved. 
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3. Improved Gluon Propagator 

Propagators and vertices necessary to calculate perturbative corrections to the staggered bilin- 
ear operators can be found in Refs. m?], ^, with the exception of the gluon propagator for the 
Symanzik-improved action. Thus we discuss only the latter here. 

The improved gluon propagator was worked out originally in Ref. [|l^]. We have found a 
convenient repackaging of the result, which we present here. This uses the notation^ 

^"=L^M' ^^:=2sin(^^/2), (3.1) 



and the following orthogonal projectors: 

<^^/iV = ^2 ' ^/iv ~ ^liv ~ ^ iiv ) (3.2) 



The inverse gluon propagator with covariant gauge fixing can then be written 



''^V — ^ •'^ p.V ~\~ S^y Cy^^y , (3.3) 



a 



^^v — S^yki^k k^ky k^ky + ^.^ , (3.4) 



where a is the gauge-fixing parameter, 

f={o)-c'P-ck^/P), (3.5) 

and 

O) = Cpi + 8crt + 8cpg , c = Crt-Cpg, c =Cpg. (3.6) 

Since we use gauge-invariant operators the matching factors are independent of a and we chose 
a = I. For the tree-level improved Symanzik action ft) = l,c = — 1/12 and c' = 0. 
Inverting S^^y , we find the improved gluon propagator to be 



[P{P -CXi)+C^X2]5j^y + c{k^ -CXl)^^v + C^{^^)^V 
f{P[k2ik2-cXl)+C^X2\-C^X,} ' 



where c = c/f, and 

xi = Tr(^) = -k^ = 2'£ k^kl , (3.8a) 

li<v 

X2 = ^ [Tr2(^) -Tr(^2)j ^ p _ (3/2)p^V (1/2)(P)^] = 3P £ P^PyPp , (3.8b) 

^ ll<V<p 

;C3 = 1 [Tr3(^) - 3Tr(^)Tr(^2) ^ 2Tr{^^)] 

= [(P)4 + 3(^)2 _ 6P(P)2 ^ g^p _ ^ 4{PfP^PSp4 . (3.8c) 



Be careful that k ^ {k ) 'his notation 
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4. Renormalization of Bilinear Operators 

The Feynman diagrams relevant to the one-loop renormalization of the bilinear operators are 
shown in Fig. |I|. Analytic expressions for these diagrams using the unimproved gluon propagator 
can be found in Refs. [Q, |8|]. It turns out that the generalization needed when using the improved 
gluon propagator is relatively simple. One simply replaces the composite gluon propagator from a 
smeared link in the /x'th direction to a smeared link in the v'th direction (which is a building block 
of the calculation) as follows: 



X a/3 



(4.1) 



The key simplification arises from the fact that this composite propagator contains off-diagonal 
terms even when the gluon propagator is diagonal (as it is for the Wilson gauge action in Feynman 
gauge), so that no new types of contribution arise when moving to the (off-diagonal) improved 
gluons propagator. One must simply evaluate the loop integrals (which is done numerically) using 
the more complicated composite gluon propagator. 




(a) 



(b) 



(c) 



S » F 




(d) (e) 
Figure 1: Feynman diagrams for bilinear operators 

Two independent calculations have been done as a cross-check. In the end, we obtain the 
one-loop renormalization factors of the bilinear operators, 



O 



Latt,(l) 



^'■^■+3(471)2 



'Yij\og{aX)+C\ 



•Latt 
U 



^Latt,(0)^ 



(4.2) 



where the superscripts indicate the order in perturbation theory, / and j are indices which run over 
the 256 bilinear operators, A is the "gluon mass" used to regularize the infrared divergences, and 
Yij is the anomalous dimension matrix. The quantities of interest are C'jf", the finite parts of the 
corrections. These have, in general, non-zero off-diagonal elements. A similar expression holds 



in continuum regularizations, and by combining with the result ( [4.2[ ) one can obtain the desired 
one-loop matching factors. Full details will be presented in Ref. Here we simply quote, in 
Tables [T] and ^, results for some representative diagonal and off-diagonal elements of C^"". 
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As noted in Ref. one can apply the mean-field improvement method of Ref. [ p2| ] to our 
HYP-smeared operators. This should reduce the size of fluctuations in the gauge links. This results 
in shifts in the values of the diagonal elements C^'^ [Q, ^ |6|]. Mean-field improved values of C^" 
are shown in the last four columns of Table |l[ 



Operator 


(a) 


(b) 


(c) 


(d) 


(ay 


(by 


icy 


idy 


(1(»1) 


41.73 


2.60 


32.82 


1.91 


31.86 


1.54 


25.59 


1.18 


(75 ® 1) 


-35.86 


-7.17 


-27.24 


-6.01 


-6.25 


-4.01 


-5.55 


-3.84 


(7m<»1) 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


(7Ai«"?5) 


-22.51 


-3.97 


-17.04 


-3.10 


-2.77 


-1.86 


-2.58 


-1.65 


(7mv«)1) 


-10.97 


-1.84 


-8.38 


-1.34 


-1.10 


-0.79 


-1.15 


-0.62 


(7mv ® ^pa) 


-34.05 


-5.19 


-25.64 


-3.82 


-4.44 


-2.03 


-3.95 


-1.65 



Table 1 : Diagonal part of representative diagonal coefficients C?/^". Note that /i , v, p , and a are all different. 
Results are given for the following four choices of actions: (a) thin links (in the fermion action and bilinears) 
with the Wilson plaquette action; (b) HYP-smeared links with the Wilson plaquette action; (c) thin links 
with the improved gluon action; (d) HYP-smeared links with the improved gluon action. The prime in the 
labels indicates that mean-field improvement is applied. Results are accurate to the number of digits quoted. 



Name 


Operator-/ 


Operator-j 


(a) 


(b) 


(c) 


id) 


CVVM 




(7m®^m) 


-3.042 


-0.351 


-2.495 


-0.321 


CVAM 


(7m ® ^ii5) 


(7m®^v5) 


0.647 


0.257 


0.609 


0.244 


CVTM 


(7m ® ^iiyi) 


(7m«'^Pv5) 


1.486 


0.280 


1.292 


0.266 


ctam 


(7mv®'?M5) 


{y^v^^ps) 


0.676 


-0.006 


0.547 


-0.003 



Table 2: Off-diagonal coefficients C'"'. The notation is the same as in Table [l|. Results are accurate to the 
number of digits quoted. 

5. Discussion 

We can use the results to compare the reduction in the size of one-loop matching factors 
achieved by different improvement schemes. For vector currents and the off-diagonal coefficients, 
the anomalous dimensions vanish and the C-"' give a direct measure of the size of the corrections. 
Looking at the results for (7^ (8) c^j), a vector current containing 3 links, one sees that the corrections 
in column (a) for unimproved fermions and glue are reduced by all choices of improvement scheme. 
Comparing the impact of applying each improvement alone, the greatest reduction is achieved by 
mean-field improvement, with HYP-smearing following close behind, but with the improved gauge 
action alone leading to a much smaller reduction. Combining all three improvements leads to the 
smallest coefficients. Similar results holds for the off-diagonal coefficients of Table ^, except that 
these are unaffected by mean-field improvement. 

For the scalar and tensor operators, which have non-zero anomalous dimensions, the coeffi- 
cients Clj' depend on the choice of infrared regularization. Thus one should consider the difference 
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between the coefficients for fixed spin and differing tastes, i.e. the differences between the first 
and second rows and between the fifth and sixth rows of Table [l[ For these differences one finds 
a similar pattern of improvements to those noted above, except that for the scalar bilinears HYP 
smearing leads to a significantly greater reduction than mean-field improvement. 

Overall, we conclude that improving the gluon action reduces the one-loop corrections some- 
what but is not nearly as effective in this regard as HYP-smearing. 
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